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A transform method for determining the flow 
generated by the singularities of Stokes flow in a 
two-dimensional channel is presented. The analysis is 
based on a general approach to biharmonic boundary 
value problems in a simply connected polygon 
formulated by Crowdy & Fokas in this journal. The 
method differs from a traditional Fourier transform 
approach in entailing a simultaneous spectral analysis 
in the independent variables both along and across 
the channel. As an example application, we find the 
evolution equations for a circular treadmilling 
microswimmer in the channel correct to third order in 
the swimmer radius. Significantly, the new transform 
method is extendible to the analysis of Stokes 
flows in more complicated polygonal microchannel 
geometries. 
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1. Introduction 

The study of low Reynolds number flows in confined 
geometries bounded by no-slip walls is of perennial 
interest in fluid dynamics [1]. It has received a resurgence 
of attention in recent years owing to new applications 
in emerging technologies in micro- and nanofluidics, 
biological fluid d3mamics and colloid science where 
the small scales render the role of inertia negligible. 
The relevant boundary value problems are linear but 
finding solutions even for the fundamental singularities 
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of Stokes flows in just two dimensions can be challenging if the geometries are complicated. 
The monograph by Pozrikidis [1] documents solutions for the fundamental Stokes singularities 
in various geometries and surveys the variety of mathematical techniques available for 
their analysis. 

This paper focuses on Stokes flow in the relatively simple geometry of a two-dimensional 
channel, but it offers a new mathematical perspective. The geometry is an important one in 
applications. Davis [2] has presented the solutions to various problems involving distributions 
of point forces, or stokeslets, in a two-dimensional channel with a view to understanding the 
blocking properties of periodic arrays of wall-attached barriers. His analysis makes use of Fourier 
transforms with respect to the variable along the channel length. Liron & Mochon [3] have found 
the Stokes flow field due to a three-dimensional force singularity at an arbitrary location and 
arbitrary orientation between two parallel plates. Section 3.5 of the monograph by Pozrikidis [1] 
describes a standard Fourier transform analysis of a two-dimensional point force in a channel. 

Although it is not the principal focus of this article, the present study was motivated 
by problems of computing the swimming dynamics of low Reynolds number organisms in 
microchannels. In a recent paper, the authors [4] have given a general prescription for computing 
the dynamical system governing the motion of a model two-dimensional circular treadmilling 
swimmer in an arbitrary domain confined by no-slip walls. By an asymptotic expansion, the 
governing ordinary differential equations are determined correct to third order in the swimmer 
radius: for a given fluid region, rather than solving a boundary value problem for the full finite- 
area swimmer, one solves boundary value problems for certain point singularities of Stokes 
flow — namely, a torque-free stresslet and an irrotational (source) quadrupole — in that domain. 
Even though this is usually simpler than solving the full boundary value problem for the finite- 
area swimmer, it remains mathematically challenging for most geometries and this challenge has 
led to the results herein. 

Our principal aim here is to present a transform approach to finding analytical solutions for 
those singular Stokes flows in the channel geometry. We adapt a new general transform approach 
to biharmonic boundary value problems in a simply connected polygon formulated by Crowdy & 
Fokas [5], who originally showcased their method for canonical semi-strip problems in plane 
elastostatics. The method differs from the usual Fourier transform approach where transforms 
are taken only with respect to the independent variable along the channel length [1,2]; instead, 
a simultaneous spectral analysis with respect to the two independent variables both along and 
across the channel is performed. It offers various mathematical advantages, not least of which is 
the fact that the approach can be generalized to any simply connected polygonal geometry. This 
work appears to be the first to apply the general method of [5] to problems of slow viscous Stokes 
flows, where biharmonic boundary value problems also arise. It showcases the method in one 
of the simplest non-trivial situations (the channel geometry) and therefore constitutes a valuable 
case study. The simplicity of the geometry also allows an instructive comparison with traditional 
Fourier transform methods detailed in an appendix. As an illustrative example, in §5, we apply 
the method to our motivating problem of treadmilling microswimmers in a channel. 



2. Stokes flows in two dimensions 

Consider a region of incompressible fluid of viscosity /x evolving according to the Stokes 
equations. It is well known that a solution of the Stokes equations for the incompressible velocity 
field (u, v) can be written in terms of a streamf unction if{x,y) with 

u= — and i; = . (2.1) 

dy dx ^ ^ 

The streamfunction then satisfies the biharmonic equation 



(2.2) 



where is the two-dimensional Laplacian. Crowdy & Fokas [5] present their approach to 
biharmonic boundary value problems in a polygon in terms of x// and its various partial 
derivatives, but there is an alternative equivalent presentation of the same method in terms of 
the so-called Goursat functions [6]. With z = x-\-iy the Goursat functions, here denoted by/(z) and 
g(z), are analytic functions of z in the fluid region with the property that the solution of (2.2) is 
represented in the form 

fiz,z) = hn[2f{z)+g{z)]. (2.3) 



The Goursat functions are familiar to fluid dynamicists using complex variable methods for 
Stokes flows [6] so we adopt this formulation and seek integral representations of the required 
physical quantities u, v, fluid pressure j) and fluid vorticity co by making use of the following 
relations of the latter quantities to the Goursat functions [6]: 



and 



4/(z) =^-ico 
-J(z) + -zf{z)+g'{z) = u-iv. 



(2.4) 



3. Point stresslet in a channel 

One important singularity of Stokes flow is the stresslet [1]. Consider a channel occupied by fluid 
at zero Reynolds number in the region — oo<x<oo,0<i/</z.A torque-free stresslet of strength 
gia ^ point zo requires that/(z) and g'{z) have the local expansions [7] 



and 



/(z) = 
/(z) = 



eVe- 



2ia 



z-zo 



+ locally analytic function 



-7^ -t- locally analytic function. 

{z - zo^ 



(3.1) 



Our objective is to find the flow induced in the channel by such a singularity. The velocity as |x| ^ 
00 is taken to vanish. We focus on this stresslet for two reasons: first, it is relevant to the study of 
swimming organisms considered later; second, we have been unable to find the solution reported 
elsewhere (although it can in principle be derived as an appropriate parametric derivative, with 
respect to the singularity position, of the solution for a point stokeslet in a channel given by 
Davis [2]). It should be pointed out that the quantity eVe^^ is just a complex constant but its 
designation has been chosen so as to be consistent with our later application of the results to a 
model treadmilling swimmer. 

It is convenient to introduce the meromorphic functions 



and 



/.(z) = .Ve-(iL)coth[|^(z-zo)] 

^^(2) = eV [e2-2^ (iL)' cosech^ [^(z - Zq)] 



(3.2) 



The functions fs{z) and gg{z) have precisely the required local form (3.1) at zq and no other 
singularities in the channel. Both functions are periodic with period 2 hi, i.e. 



fs(z + 2hi)=fs(z) and g',{z + 2hi)=g',(z). 



(3.3) 



The derivative /g(z) is also 2 /zi-periodic and this will be useful later. 



We now pose the ansatz 

/(2)=/,(z)+/r(z) and g'{z)=g',{z)+g'j,{z), 



(3.4) 



where /r(z),^^(z) are analytic in the channel and vanish as |x| ^ oo; note that the second term in 
the definition of gg(z) has been chosen to ensure that 



u — iv = —f{z) + zf{z) + g\z) 0, as |x| oo. 



(3.5) 



It should be emphasized that it is not necessary to introduce the hyperbolic functions (3.2) and 
one can proceed instead by defining /s(z) and gg{z) to have the simple rational function forms 
in (3.1). However, we have found that the periodicity properties (3.3) greatly simplify the algebraic 
manipulations to follow. 

The no-slip conditions on the two channel walls are 



— /(z) + zf(z) + g\z) = 0, on z = z and z = z — 2/zi. 
On substitution of the ansatz (3.4), we find that, on z = z, 

-W) + zfR(z) + g'R(z) =m - z/;(z) - g'siz), 

and, onz = z — 2hi, 

-J^) + (z - 2hi)fj,{z) + g'j,{z) =Js{z - 2hi) - (z - 2M)/;(2) - ^^(z). 



(3.6) 



(3.7) 



(3.8) 



where known functions appear on the right-hand sides of (3.7) and (3.8). The Schwarz conjugate 
of an analytic function /z(z) is defined by 



h(z) = h(z). 



(3.9) 



Following the general formulation in [5] the two functions /r(z) and ^^(z), which are analytic in 
the channel and decaying in the far-field, will be represented by 



where pi{k) and p2{k) are defined by 



P2{k)e'^^dk 



]■ 



Pi(fc)= 



/R(z)e-*Mz and P2{k) = 



— oo+i/z 
oo+i/z 



/R(z)e-'*Mz, 



and, similarly. 



1 

2^ 



P3{k)^^^dk^ 



P4{k)^^^dk 



where psik) and p^ik) are defined by 



P3{k) = 



g^(z)e ^'^^dz and P4{k) = 



-oo+i/z 



oo+i/z 



^Uz)e-'''dz. 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



The set {pn{k)\n = 1, 2, 3, 4} is called the set of spectral functions [5,8], and there are relations between 
these spectral functions referred to as global relations. In this case, there are two global relations 
given by 

P2{k) = -pi(k) and P4{k) = -p3{k), keR. (3.14) 

In this case, these relations can be viewed as a consequence of Cauchy's theorem: for example, 
the condition p2{k) = —pi{k) is equivalent to the statement that the integral of/R(z)e~^^^ around 
the channel is zero because, for /c e R, the latter function is analytic inside the channel. On use of 



these relations (3.10) can be written as 

1 



fR{z) = 



In 



0 



or, on making the substitution k\-> —kin the second integral, 

1 poo |- -. 

/r(z) = ^ [Plik) e*^ + pi(-fc) e-*^] dk. 

Similarly, we can write 

1 r 1 

gR{z) = ^ [P3{k) + P3(-^) e-^^^] dk. 



(3.15) 



(3.16) 



(3.17) 



It is worth emphasizing that it is this spectral representation of the analytic functions /^(z) and 
g^(z) that can be extended to arbitrary simply connected polygons [5,8,9]. 

(a) Preliminary observations 

The representations of /j^(z) and ^^(z) above involve only real values of k. Hence, for real k, 
complex conjugation implies 



Pi(k)-- 



/R(2)e"^Mz. 



If we \etk\-> —k we find 



fR(z)e-'^'dz = pi{-k). 



(3.18) 



(3.19) 



The integral on the left-hand side will appear in the transform of the boundary conditions on the 
wall at y = 0. 

Similarly, on writing x = z — ihin the integral expression for p2{k), we can write 



/R(x + i/z)e-^^(^+^^) dx. 



P2{k) = - 

On taking a complex conjugate for real k, 

P2{k) = - fR{x + ih) e^^(^-^^) dx = -e^^ 



/K(x + i/z)e^^Mx, 



Now making the substitution x = z — ihin the last integral produces 

f oo+i/z 
-oo+i/z 



P2{k) = -e- 



2kh 



fR{z)e'''dz. 



Hence, mapping k\-> —k,it follows that 

rOO+i/z 



Mz)e-'^'dz = -e'''p2(-k) 



2kh- 



-oo+i/z 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



Again, the integral on the left-hand side will appear in the transform of the boundary conditions 
on the wall aty = h. 



(b) Relations between spectral functions 

The boundary conditions on the channel walls allow us to deduce more information about the 
spectral functions. To proceed, we multiply the boundary condition (3.7) by e~^^^ and integrate 



along the real z-axis: 



where 



/R(z)e-*Mz + 



z/^(z)e-'*^=^dz + 



;(z)e-"^Mz = Ri(fc), 



[f,(z)-2/;(z)-^;(2)]e-'*^Mz, 



The function Ri (k) is known. On use of the result, available from an integration by parts. 



z/^(z)e-"=^dz = i-[ifcpi(fc)] = 



9[fcpi(fc)] 
dk ' 



together with (3.19), we find 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



Similarly, if we multiply the second boundary condition (3.8) by e and integrate over the upper 
boundary, we find 



e^^^Pii-k) - iir[^kp2{k)] + {2ih)[ikP2{k)] - p^ik) = Ri{kl 



dk 



where 



Mk): 



oo+i/z 



-oo+i/z 



\fs{z - Ihi) - (z - 2/zi)/;(z) - g',{z)\ e-''' dz 



(3.28) 



(3.29) 



Equations (3.27) and (3.28) can be analysed, together with the global relations (3.14), to determine 
the spectral functions. 

Indeed, on use of (3.14) equation (3.28) becomes 



- e^''^M=k) - + IhkMk) + p3(fc) = R2ik). 

On subtraction of (3.27) from (3.30), we find 

- (e^*^'' - + 2hkpi{k) = R(k), 

where 

R{k)^R2{k)-R^{k). 
After taking a complex conjugate and letting k \-> —k, we find 



(1 - e-^^^)pi{k) - lhkpi{-k) = R(-k). 



Elimination of pi{—k) from (3.31) and (3.33) produces 

, -2 hkR{k) + (e^^^ - 1)R(^ 

Pi (k) = 7^ . 

4{smh^{hk) - h^k^) 

The function R{k) can be computed explicitly using residue calculus to give 



(3.30) 

(3.31) 
(3.32) 

(3.33) 
(3.34) 



Pi{k) = - 



nieV 
2A 



2hk 



-2kh 



where 



+ ik{Zo - Z^) [2Me^(2a-fczo) _ ^-i{2a+kz-o)^^2kh _ -^^j 

A{k) = sinh^ {hk)-hV. 



(3.35) 



(3.36) 



Given pi{k), the spectral function p^{k) follows from (3.27) together with 



Riik)-- 



27ri£yei(2«-fc2o) 

p2fc/7 _ 1 



-2kh 



- 1 + 2/c/z 



1 



+ i/c(zo - Zo) 



(3.37) 



which also follows from a residue calculation. All physical quantities are now available from 
the representations, dependent only on/s(z),^g(z) and the associated spectral functions pi{k) and 
Ri (k), given by 

^-ioj = 4/(z) = 4/;(z) ^-T \ikpi(k) e'^^ - ikpi(-k) e"^^^] dk 

/X Jo L -I 



and 



Jo 

:-^+z/;(z)+^;(z) 
1 r 

27r Jo 
— i(z — z 



[Riik) 

l)k[pl(k) 



e*^ + Ri(-/c)e-*^ 



e'*^^-pi(-fc)e-*H 



+ PiW(e" 



-ifcz 



e-*^) + Pi(-fc)(e' 



ifcz 



e^^^)] di^. 



(3.38) 



This is the required solution. Significantly, these integrals are readily computed numerically 
owing to the exponential decay of all integrands as /c ^ oo. As /c ^ 0, 

b 



where 



and 



671 



eV[h + i(zo — Zo)] sin 2a 



c = - ^eV{2h^ cos 2a + 3[h + i(zo ■ 



■zo)][2Im(zo e 



2ia\ 



/ze-2^^]}. 



(3.39) 



(3.40) 



Note that h is real-valued and c becomes purely imaginary if yQ = h/2 or sin2a = 0. 
Expressions (3.40) can be used to verify that all integrands in (3.38) are regular at k = 0. 
Formula (3.39) is also useful in facilitating explicit computation of the limiting behaviour of the 
integrands in any numerical evaluation procedure. Finally, it should be pointed out that although 
the integrals defining all physical quantities in (3.38) are regular at k = 0 the integrand in the 
expression (3.16) for/]^(z) is singular at A: = 0 but this is rendered immaterial by a freedom to 
redefine /(z) i-^/(z) -\-pz-\- q and g\z) \-> g\z) + q for constants p eR,q eC without affecting any 
of the physical quantities. 

It is instructive, for purposes of comparison, to solve the problem of a point stresslet in a 
channel using a more traditional approach involving a Fourier transform only with respect to x. 
This separate analysis in appendix A also affords us a numerical check on (3.38). 



4. Point quadrupole in a channel 

It is straightforward to adapt the approach to other singularities. Indeed, the solution for a stokeslet 
singularity [1] in a channel can be derived by simple adaptation of the above approach and 
compared with the result of Davis [2], who employed standard Fourier transform methods akin 
to those used in appendix A for the stresslet problem. For use in §6, consider the separate problem 
in which an irrotational quadrupole is situated at zo with the corresponding function/(z) analytic 
at Zo and g\z) behaving like 

2eVe^^'' 

g (z) = -r^ + locally analytic function. (4.1) 

A derivation of the flow generated by this singularity closely follows that of §3: the functional 
forms in (3.38) are exactly the same, the only difference is that/s(z),gg(z) must be adapted thereby 
producing correspondingly different spectral functions pi (k) and Ri (k) (it is convenient here to 



use the same notation as in §3). But it turns out that we can bypass a full analysis by noticing that 
the functions /(z) and g\z) relevant for the point quadrupole problem can be derived from those 
of §3 by computing the mixed parametric derivatives 



/(z)^^ and g'iz)^'"^'^'^ 



dZodZo 



dZodZQ 



(4.2) 



This generation of the flow due to a quadrupole by taking parametric derivatives with respect 
to the singularity position is equivalent to use of the Vq operator by Davis & Crowdy [4] in 
their analysis of low Re3molds number swimmers in confined domains (the application to be 
considered in §5). It follows that, for the point quadrupole problem. 



/s(z) = 0, g^^{z) = 2eVe''-[^] cosech^ 



7t{z-Zo)' 

2h 



coth 



Tziz-ZoY 

2h 



(4.3) 



with spectral functions in the integral representations for the corresponding /j^(z) and ^^(z) now 
given by 



and 



Ri(fc) = 



i 

2A 



Q2kh _ I 



(4.4) 



5. Application: microswimmers in a channel 

One application of the foregoing analysis is to modelling the low Re3molds number swimming 
of micro-organisms in a channel. In 1963, Lord Rothschild [10] conducted an experiment in 
which a sample of bull sperm in a very viscous fluid was placed between two glass plates; after 
time he measured the distribution of the sperm between the two plates and found a tendency 
for them to aggregate at the channel walls. The experiment was repeated recently by Berke 
et al. [11] and a hydrodynamic explanation of the wall attraction phenomenon proposed. But 
more general dynamical behaviour of low Reynolds number swimmers near walls has been 
observed, including steady translation and periodic motions along the walls [12,13]. To rationalize 
this behaviour using simple mathematical models, Crowdy & Or [7] have proposed the study 
of a simple circular 'treadmiller' comprising a cylindrical circular swimmer with an imposed 
tangential velocity profile; this profile actuates motion in the spirit of an 'envelope model' of 
surface ciliatory motion [14]. Crowdy & Or [7] showed that, in free space, such a swimmer 
has the singularity distribution comprising a torque-free stresslet and a superposed irrotational 
quadrupole; as an approximation, they studied the motion of such a singularity combination 
near a straight no-slip wall and found good qualitative agreement with numerical [12] and 
laboratory [15] experiments. Davis & Crowdy [4] have since found the general evolution 
equations for the same swimmer correct toO(e^), where e is the swimmer radius. We now use the 
latter equations to study the dynamics of such a treadmilling swimmer in a channel. 

The swimmer, which is oriented at some angle a{t) to the real axis, has an imposed tangential 
surface velocity of the functional form 2 sin 2(0 — a{t)), where 0 is an angular coordinate around 
the swimmer boundary. The interaction of the Stokes flow it generates with any nearby solid 
walls generally results in dynamical evolution of its centre Zo{t). A schematic of such a swimmer 
located in a channel of width h is shown in figure 1. 

The evolution equations given in Davis & Crowdy [4] have the following form. Suppose 
\I/R denotes the streamfunction of the reflected velocity generated by the wall interaction of the 
stresslet considered in §3, where 



fR = lm[zfR{z)-\-gR{z)l 



Figure 1. Circular treadmilling swimmer, of the type proposed in [7], centred atzo(f) at orientation a(t) in a channel of width 
h. The radius of the swimmer is € < 1. Its evolution equations, correct to 0(€^), can be determined from the results in [4] 
using the transform results of §§3 and 4. (Online version in colour.) 



with 



and 



fR{z) =fo{zo,a) +/i(zo,a)(z - zo) +/2(zo,a)(z - Zof + • 
gR{z) = goizo^oi) -\- gl{zo,a){z - Zo) + • • • 



(5.2) 



Suppose too that the streamf unction ^jfR for the reflected velocity generated by the point 
quadrupole of §4 can similarly be determined in the form 



xlfR = lm[zfR{z)-\-gR{z)l 



where, near zq. 



fR{z) =fo{zo,a) +/i(zo,a)(z - zq) +/2(zo,a)(z - Zof + • 
and g'j^iz) = go{zo,a) gi{zo,a){z - zq) + • • • 

The evolution of the swimmer, correct to 0{e^), is given [4] by the solution of the system 
dzoit) 



(5.3) 



(5.4) 



and 

where 



dt 
da 
'dt 



= e(us + ivs) + 6^(WQ + wq) + €^(UA + iva) 



(5.5) 



Us + ivs = -fo{zo,a) -\- Zofi{zo,a) go{zo,a), 



and 



+ i^Q = -fo{zo,a) + zofi(zo,a)-\rgo{zo,a), 
ua + iVA 
Qa 



^2s = -2Im[/-i(zo,a)], 
QQ = -2\vcv\fi{zQ,a)l 



2f2{zo.a\ 

-Im[2/i(zo, -7r/4)]Im[2z^/2(zo,a) +gi(zo,a)] 
- Im[2/i(zo,0)]Re[2z^/2(zo,a) +^i(zo,a)]. 



(5.6) 



The main point here is that the analysis of §§3 and 4 allow all the quantities (5.6) appearing 
in the approximate evolution equations (5.5) to be calculated and, thus, the swimmer motion 
determined. 

Figure 2 shows some typical swimmer trajectories for a{0) = 7r/4 and x(0) = 0,i/(0) = 0.1(0.1)0.9 
for h = l. Consistent with the experimental observations of wall aggregation [10,11], it is found 
that the swimmer is attracted to the nearest wall of the channel with only swimmers exactly 
positioned on the channel centre-line remaining stationary. A range of more general initial 
orientations were explored and all swimmers off the centre-line are similarly found to be attracted 
to the nearest wall. We found no initial conditions for which the swimmer propagates along the 
channel. It is important to recall that these evolution equations are approximate and only correct 
to 0(e^); it would be of interest, for purposes of comparison, to compute the full unapproximated 
evolution equations but this would probably require a numerical treatment using, for example, 
boundary integral methods [1]. 




-1.0 -0.5 0 0.5 1.0 

Figure 2. Swimmer attraction to channel walls with /? = 1: trajectories for € = 0.03 with initial conditions a(0) = 7r/4for 
x(0) = 0,y(0) = 0.1(0.1)0.9. The swimmer remains stationary when y(0) = 0.5. The evolution has been computed correct to 
0{€^) from the asymptotic equations given in [4]. (Online version in colour.) 



6. Discussion 

Although our motivating application was the dynamics of low Reynolds number swimming 
organisms, the principal contribution of this paper is mathematical: a novel transform approach 
to finding the flows due to basic point singularities of the Stokes equations in a simple two- 
dimensional channel has been given. Other types of Stokes flow in the channel can be studied 
using adaptations of the same approach. 

The solutions can be derived using alternative approaches, such as traditional Fourier 
methods, but the method here, based on a general formulation by Crowdy & Fokas [5], is 
significant in that it can be generalized to finding flows in more complicated simply connected 
polygonal domains (a similar solution scheme for the Laplace equation in a polygon has been 
given by Fokas & Kapaev [9]). In those cases, the determination of /(z) and g\z) will follow 
similar steps: a convenient form of the singular terms fs{z) and gg{z) must be chosen based 
on considerations from the geometry being analysed and the integral representations for the 
correction terms /]^(z) and ^^(z) will be similar except that every additional side of the polygonal 
domain will introduce a new spectral function. The set of spectral functions will again satisfy 
a global relation which must be analysed to determine any unknown spectral functions [5,8]. 
Here, the global relations are manifested in equations (3.14). In many cases, the analysis of 
these involves the solution of a Riemann-Hilbert problem for the spectral functions [5,8] but the 
problem in a channel geometry studied here does not require any such Riemann-Hilbert analysis 
and the spectral functions can be determined using only algebraic manipulations. Other boundary 
value problems for which a Riemann-Hilbert analysis can be avoided, and only algebraic 
manipulations used, have been reported by Fokas [8]. The latter monograph also provides an 
interesting survey of how the transform method for Stokes flow used here fits into a unified 
framework for boundary value problems developed in recent years by Fokas and collaborators. 

Extensions of our approach to finding Stokes flow solutions, and studying microswimmer 
dynamics, in more complicated polygonal channel geometries are currently under investigation. 

Acknowledgements. D.G.C. acknowledges financial support from a research grant from The Leverhulme Trust in 
the UK. 



Appendix A. Standard Fourier transform approach 

In this appendix, the problem for a point stresslet in a channel is solved using standard Fourier 
transform methods. Suppose a stresslet of strength e^^'^ is at zq = iV. The streamf unction x// can be 



written as 



\l/{x,y) = Im 



2ia (Z-^o 



z-zo 
= sin 2a — 2 sin 2a 



x2 + (y-y)2 



-2 cos 2a 



Lx2 + 



iy-Y) 



(y-v)2 



+f, (A 1) 



odd in X 



where is regular in the channel. It is convenient to introduce image singularities in order to 
satisfy the no normal flow conditions on the walls. Now 



and 

oo 

i: 

m=— oo 

We therefore write 



(y - 2mh - Yf (y - 2mh + Yf 

- + (y - 2mh - Yf- + (y - 2mh + Y)^ 



y - 2mh -Y y - 2mh + Y 



+ (y - 2m/z - y)2 x^ + {y- 2mh + y)^ 



= 0, ony = 0,/z (A2) 



= 0, ony = 0,/2. (A3) 



where 



and 



Then 



Now 



V^(x,y) = sin2a + 2sin2a if^^\x,y) -2cos2a if^^\x,y), 



f^'\x,y) = - E 

m= 
oo 

= E 

m=— c 



(y - 2mh - Yf (y - 2mh + y)^ 



x2 + (y - 2m/z - y)2 x2 _^ _ 2;tz/| + y)2 
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x2 + (y - 2m/z - y)2 x2 _^ _ 2/^/| + y)2 
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f^') = 0 = ^^^), ony = 0,/z. 



(0) 



_x2 + (y _ 2m/z - y)2 x2 _^ _ 2;tz/2 + y)2 _ 



^ poo 

= V e"^^[cos % - 2m/z - y) - cos k{y - 2mh + y)] dk 
Jo 



= 2 E 

m=—oo 
oo 

= 2 E 

m=— OO 
oo 

= 4- E 



e ^ sin /c(y - 2m/z) sin kY dk 



e sin /cy cos 2A:m/2 sin /cyd/c 



e-^^ sin % 5(2/:/z - 2n7r) sin /cY d/: 
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e— /^sin(!!|^)sin(!f), 



(A 4) 



(A 5) 



(A 6) 



(A 7) 



(A 8) 
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3> 



by use of the periodic <5 -function. Similarly, 
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In 
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Hence 
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nnY 



if 
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/n7Ty\ ( n7iY\ 



(0) 



The boundary conditions on the channel walls now take the form 



and 



ay 



ay 



2 oo 



nnY 



Next, we introduce the Fourier sine and cosine transforms defined by 



Conditions (A 7) imply that 
The identity 

implies the results 



f cos Xxdx and ^ -- 

= 0 = ip^^\ ony = 0,/z. 
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which can be used to take the transforms of (A 11) 
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(A 9) 



(A 10) 



(All) 



(A 12) 



(A 13) 



(A 14) 



(A 15) 



(A 16) 



The Fourier transform of the biharmonic equation satisfied by the streamfunction implies that 



dy2 



a2| 0^ = 0. 



(A 17) 



Independent solutions of this ordinary differential equation that vanish on the channel walls at 
y = Q,h are 



/ y\ sinhAy ^ y / cosh Ay sinhAy\ 
V h) smhXh /z IcoshA/z smhXh)' 



(A 18) 



However, inspection of (A 16) shows that, for each function, odd /even values of n yield the 
symmetric /antisymmetric parts with respect to the midplane y = h/2. Hence we let 

(h — y) sinh Ay + y sinh[A(/z — y)] 



and 
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(A 19) 



sinh(A/z) — Xh 

where and are to be determined. Note that both functions remain 

bounded as A ^ 0. Substitution of (A 19) into (A 16) then yields 
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and 
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Finally, the inverse cosine and sine transforms give 



^W(x,y) = - 
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The velocity components can be found by differentiation of these expressions. It can be verified 
numerically that the values obtained are the same as those computed from (3.38). 
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